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Abstract. Dimensional continuation is employed to compute the energy loss rate 
for a non-relativistic particle moving through a highly ionized plasma. No restriction 
is made on the charge, mass, or speed of this particle, but it is assumed that the 
plasma is not strongly coupled in that the dimensionless plasma coupling parameter 
g = e 2 K D /4irT is small, where k d is the Debye wave number. To leading order in 
this coupling, dE/dx is of the generic form g 2 ln[<? 2 C]. The prefactor of the logarithm 
is well known. We compute the constant C under the logarithm exactly. Our result 
differs from approximations given in the literature, with differences in the range of 
about 20% for cases relevant to inertial confinement fusion experiments. 
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1. dE/dx and the Coulomb Log 

The stopping power of plasma component b for projectile p is of the form 

dE b = e 2 p e\ n b = ^ — In A (1) 

dx Air m b Vp Air f3 b m b Vp 

where the Coulomb logarithm \nA b involves a ratio of short- and long-distance length 
scales. To compute lnA^, we employ the method of dimensional continuation [1]. To 
introduce this method, we consider the Coulomb potential <fr u (r) of a point source in v 
spatial dimensions: (j> v (r) ~ l/r u ~ 2 . Clearly the long and short distance behavior of <f> v 
depends on the spatial dimensionality v. In high u, short distance (hard) interactions 
are accentuated, while in low v the large distance (soft) physics predominates. 

For Coulomb interactions, v = 3 is special in that neither hard nor soft processes are 
dominant. For v < 3, the soft physics is predominant, and for v > 3 the hard processes 
are dominant. The energy loss for v > 3, dE B /dx, is obtained from the Boltzmann (B) 
equation, and it contains a pole {y — 3)" 1 that reflects an infrared divergence in the 
scattering process when v — > 3 + . The energy loss for v < 3, dE hB /dx, is obtained from 
the Lenard-Balescu (LB) kinetic equation, and it contains a pole (3 — v)~ l that reflects 
an ultraviolet divergence when v — > 3~ . The stopping power to subleading order is 
therefore 

dE , ( dE LB dE B , 

(2) 
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and it is completely finite. Hence the two poles must cancel. The dependence of the 
residues of the poles on v brings in a logarithm of the ratio of the relevant short- and 
long-distance length scales, which is precisely the Coulomb logarithm. 

2. Collective Excitations: Lenard-Balescu Equation for v < 3 



The soft physics is described to leading order in the plasma density by 

^/a(Pa) = -£ 6 ^- J a<>> 

which is the Lenard-Balescu kinetic equation for plasma species a and b, where 
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(2tt)" (k 2 ) 2 | e (k,v a -k)| 2 y {2-KhY "V~ -'J ~ [d Pb dpa 
with v a = p a /m a and v& = p b /m b . The collective behavior of the plasma enters through 
its dielectric function 



d v p c 



1 



k~/ e (Pc), t/-0 + .(5) 



(2iih) v lo — k ■ v c + irj dp c 
The rate of kinetic energy loss of species a to species b is given by 
dS^_ = I d v p a p 2 a d 
dt ' J {2nhy 2m a dp a ah ' 1 ' 

We evaluate this for the case in which species a is a single projectile of mass m p and 
velocity v p , / a (p a ) = (2irh) u 5^ (p a — m p v p ) , and the distribution function fb(Pb) for 
plasma species b is Maxwell-Boltzmann at temperature TJ, = With dx = v p dt, 
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where M p b = m p + nib is the total mass, F(v) = k 2 [e(k, /a>) — 1 ], and = 
til^J (3brrib/2'n v exp(— ^ra^v 2 j2\ Here Q v is the area of a unit sphere in v dimensions 
and K is an arbitrary wave number whose dependence cancels in the limit (2). 

3. Hard Collisions: Boltzmann Equation for v > 3 



Hard collisions in the plasma density are correctly described by the Boltzmann equation, 
which gives 



dEl 
dx 



d u Pb , ( , f, 1 



1 2 2 



(8) 



where v p f, = |v p — V&I, and da p f, is the full quantum-mechanical differential cross section for 
scattering of the projectile (p) from the initial velocity v p = p p /m p to the final velocity 
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Vp off a plasma particle (6). Straightforward kinematical manipulations exploiting the 
axial symmetry of the scattering produce the form 



dEf 
dx 
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(2nh)» JU ^ u/ 2p 2 M pb 

in which P is the total momentum of the center of mass, p is the relative momentum 
in the center of mass, and q is the momentum transfer. 

The classical cross section in v dimensions is do~p b = Q u -2 B u ~ 2 dB, where B is the 
classical impact parameter. Some calculation gives 
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(10) 



with rripb = m p mb/Mpb being the reduced mass. Placing the result (10) in Eq. (9) yields 
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Making the decomposition / da p bq 2 
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J dOp h q 2 + j(da p b — da pb ) q" expresses 



dEf /dx = dE bc /dx + dE bQ /dx, where dEf Q /dx is the quantum correction to Eq. (11). 
The integral / (da p b—da pb ) q 2 is most easily evaluated by first calculating / (da p b—do-p b ) q 2 , 
where dap\ is the Born approximation to do~ p b, and then subtracting the contribution 
J(,do~pb — do~p b ) q 2 . Inserting the correction v p b J(da p b — do~p b ) q 2 into Eq. (9) yields 
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(12) 



where i/j is the logarithmic derivative of the gamma function, Re denotes the real part, 
at = fibiribV 2 , and r]b = e p eb/A7rhv p . 



4. Results 



The total stopping power is the sum of the contributions from large- distance collective 
excitations dE hB /dx and from short-distance hard collisions dE B /dx, that is, the sum 
over species b of Eqs. (7), (11), and (12). The poles at v = 3 and the \nK terms cancel. 
Our result for dE/dx is generically of the form n (lnn + C) in the plasma density n, and 
it is accurate to all orders in the quantum parameter rj b - Figures 1 and 2 illustrate our 
result with an example that is relevant to the DT plasmas in laser fusion capsules. 
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Figure 1. The energy E(x) (in MeV) of an a particle with initial energy Eq = 
3.54 MeV vs. the distance x (in /im) traveled through an equal molal DT plasma. 
Note that E(x) is obtained by inverting x = J Eo dE {dE / 'dx)~ l , where the stopping 
power dE/dx has been expressed as a function of energy. The plasma temperature is 
T = 3 keV and the electron number density is n e = 10 25 cm -3 . The plasma coupling 
is small, g = 0.011, and so our calculation (BPS) is essentially exact. Our result is 
shown by the solid curve. The work of Li and Petrasso [2] is often used in laser fusion 
simulations. Their result (LP) is shown by the dashed line. Note that the difference 
in the total ranges between our result and that of Li and Petrasso of about 5 \xm is a 
little larger than 20%. 
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Figure 2. The a particle dE(x)/dx (in Mev/fim) vs. x (in /jm) split into separate ion 
(spiked curves) and electron components (softly decreasing curves). The area under 
each curve gives the corresponding energy partition into electrons and ions. For our 
results (BPS), the total energy deposited into electrons is E e = 3.16 MeV and into 
ions is Ei = 0.38 MeV, while LP give Ef = 3.11 MeV and Ej F = 0.43 MeV. These 
energies sum to the initial a particle energy of Eq = 3.54 MeV. Note that BPS has a 
longer a particle range and deposits less energy into ions than LP. Both observations 
would tend to make fusion more difficult to achieve for BPS than for LP. 




